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STUDY OF HIGHLY SWEPTBACK WINGS BY THE 
FREE VORTEX SHEET METHOD 

By 

C. Subba Reddy ^ and Farhad Ghaffari^ 

ABSTRACT 

The aerodynamic characteristics of highly sweptback wings with separa- 
tion induced vortex flows have been numerically investigated using the free 
vortex sheet method, developed by Boeing Company, under a contract with 
NASA/Langley Research Center. The models studied included delta and straked 
wings, and wings with leading edge extensions. Also, PAN-AIR code has been 
used to design a fixed leading edge extension to a thick delta wing. 

The theoretical results predicted have been compared with the experi- 
mental data wherever available, and the code capabilities and limitations 
explored. The fuselage effects on the aerodynamic performance have also 
been considered in some cases. 


INTRODUCTION 

This report briefly describes the research conducted under grant NSG 
1561 during the period September 1, 1981 to August 15, 1982. In this work, 
the latest version of the free vortex sheet (FVS) method of Boeing Company 
(refs. 1 and 2) has been mainly employed to study various configurations not 
covered by the previous reports (refs. 3-6), with a view to determining the 
code capabilities and limitations. Also the recently incorporated code 
capability for modeling wings with multiple vortex systems has been utilized 
on double delta wings. The quasi-vortex lattice (QVL) method of Mehrotra 
(refs. 7 and 8) and the vortex lattice method with the suction analogy (VLM- 
SA) of NASA/Langley Research Center (refs. 9-11) which were extensively used 
in previous studies (refs. 3 and 4) have not been employed in this investi- 
gation. However, another code, PAN-AIR (ref. 12) has also been used to a 

^Assistant Professor, Department of Mechanical Engineering and Mechanics, 

Old Dominion University, Norfolk, Virginia 23508. 

^Research Assistant, Department of Mechanical Engineering and Mechanics, Old 
Dominion University, Norfolk, Virginia 23508. 



limited extent. The models studied included thin delta and straked wings 
with fuselage attached. Also included are thick, and leading-edge flapped 
wings. The details of these planforms are given in table 1. In the follow- 
ing sections, the results are discussed and some of the code capabilities 
and limitations are evaluated. 
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NOMENCLATURE 

aspect ratio 
wing span 
local wing span 
local wing chord 

mean aerodynamic chord 
wing root chord 

drag due-to-Iift coefficient 
lift coefficient 
pitching moment coefficient 
net lifting pressure coefficient 

drag 

free vortex sheet 
Mach number 
body axis coordinates 
angle of attack 


RESULTS AND DISCUSSION 


In this section, the results obtained by using the numerical codes are 
compared with the experimental data, wherever available, and the code 
capabilities and limitations are discussed. A summary of the configurations 
investigated is presented in table 1. The angle of attack range over which 
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the code is employed, the type of method used, and whether convergence is 
obtained or not, are also indicated in the table. 

Figure 1 shows the delta and double delta wing models which are inves- 
tigated using the FVS method (ref. 1). These are 63.5* delta and 78.7*/ 
63.5* double delta wings. The aerodynamic characteristics predicted by the 
method for the delta wing-body combination are compared with the experiment- 
al data (ref. 13) for various Mach numbers in figures 2-4. In figure 4, the 
results obtained for the wing only are also shown. As the figures indicate, 
the agreement between the predicted results and the data is not very good. 

In figures 5 and 6, the spanwise pressure distributions at two chordwise 
stations, obtained from different types of modeling of the delta wing by the 
FVS method, are compared with the experimental data. The data is not shown 
for the lower surface as it is not available. The predicted results and the 
data near the apex agree more favorably than they do near the trailing edge. 
In both cases, the predicted pressure peAk is much higher than the actual 
one. 

It is suggested in Boeing's instructional manual (ref. 2) that a carry- 
over lifting system, which extends the wing lifting system into the body, 
has to be used in order to model the body effects. However, when no carry 
over lifting system is used, the program still works and provides the re- 
sults, which seem reasonable, for this wing as shown in figures 5 and 6. 

When the carry-over lifting system is used, it coincides with the upper body 
surface because the body exists only on the underside of the wing. The 
pressure values obtained from this model are unrealistic on and near the 
body as can be seen from the figures. Therefore, this type of modeling has 
not been used in further investigations. However, when the carry-over lift- 
ing system and the upper body surface are separated by a small distance 
(0.01 compared to = 250) by lifting up the upper body surface, the 
results are dramatically improved. The results obtained by this modeling, 
and the modeling without any carry-over system, are essentially the same as 
shown by the figures. However, the solution convergence is faster and, 
hence, less computational time is needed in the case of modeling with no 
carry-over lifting system. Therefore, throughout this work, this type of 
modeling has been used whenever the body effects are to be included. 
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Figures 7-10 show the spanwise pressure distributions at various chord- 
wise stations, angles of attack and Mach numbers. Experimental data is 
compared with the predicted results, wherever available, and found that 
there is not a good agreement between them in the outboard region of the 
wing. 

The 78.7“/63.5* double delta wing is modeled using two separate vortex 
systems on inboard and outboard leading edges. The predicted longitudinal 
aerodynamic characteristics compare favorably with the experimental results, 
as indicated by figure 11. Spanwise pressure distributions at different 
locations for various angles of attack are illustrated in figures 12 to 23. 
The predicted pressure distribution on the forward portion of the wing, 
before the break in the leading edge occurs, agrees reasonably well with the 
experimental data; in the aft portion, the agreement is comparatively not 
good. However, the method predicts the two pressure peaks as the data 
shows . 

A thick, round-edged delta wing glider, also called DM-1 (ref. 14), was 
originally studied in 1946 with the aim of developing a supersonic airplane. 
The DM-1 glider has approximately a delta planform and NACA-0015-64 airfoil 
sections with an aspect ratio of 1.8 and a 60® swept-back leading edge. 

The flat DM-1 

glider with Flat Leading Edge Extension (FLEE) and vertical fin removed, is 
modeled by the FVS method (ref. 1) as shown in figure 1 ^. The aero-dynamic 
characteristics obtained are compared with the data at a = 15® and M = 0.0 
in table 2. Further investigation of DM-1 using the FVS method is in 
progress. 

Another thick wing with camber and twist is being investigated using 
PAN-AIR CODE (ref. 12) for design purposes. The objective of this study is 
to develop the technology by which a Fixed Leading Edge Extension (FLEE) can 
be added to a wing without affecting its performance at cruise. The FLEE 
device, to be designed, would lie along the stagnation stream surface for 
the wing at its cruise angle of attack. At angles of attack greater than 
for the cruise, vertical flows would be generated by the device and based on 
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the results for the DM-1 glider the increase in drag can be minimized. For 
this purpose the flow fields at different spanwise stations are analyzed and 
stagnation points and surfaces determined. Such analysis is shown in figure 
for a particular airfoil section at the wing root and for a designed 
angle of attack of 6.16®. Further investigation of this study is in 
progress. 

CONCLUSIONS 

The aerodynamics performance of low-aspect ratio sweptback wings with 
vortex flows has been investigated using a numerical code, developed by 
Boeing Company, and sponsored by NASA/Langley Research Center. Wings of 
different planforms have been studied and the predicted results compared to 
the existing data, wherever possible, in order to determine the code 
capabilities and limitations. Another Boeing Company code, PAN-AIR, has 
also been used to develop a fixed leading edge extension to a thick delta 
wing. 

The longitudinal aerodynamic characteristics predicted by the free 
vortex sheet method for the delta wing-body combination do not compare 
favorably with the data. However, there is a fairly good agreement in the 
case of a double delta wing with no fuselage effects considered. Also the 
method provides results for the above delta wing with no body effects 
considered which agree better with the data than those when the body is 
included . 

In the case of spanwise pressure distributions, the predicted results 
and the data agree more favorably near the apex than they do near the 
trailing edge. However, in both cases, the predicted pressure peak values 
are higher than the actual ones. The method with double vortex modeling 
predicts two pressure peaks as the data shows in the case of double delta 
wing . 
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Table 1. Wing 


Wing Description Modeling Details 

63,5** flat delta 

wing, A = 2.0 Without near wake 

63,5** flat delta Without near wake or 

wing-body com- carry-over lifting 

bination, A = 2.0 system 

63,5® flat delta Without near wake or 
wing-body com- carry-over lifting 

bination, A « 2.0 system 

63.5® flat delta With near wake, and 
wing-body com- without carry-over 

bination, A = 2.0 lifting system 

63.5® flat delta Without near wake, and 
wing-body com- with carry-over lifting 

bination, A = 2.0 system coinciding with 

upper body surface 

63.5® flat delta Without near wake, and 
wing-body com- with carry-over lifting 

bination, A = 2.0 system and the upper 

body surface separated 
by a small distance 
(0.01); discontinuity 
exists between body and 
wing upper surfaces 


00 


ions modeled by the FVS method. 


Angle of 
Attack 

Mach 

Number 

Method of Solution 

Solution 

Converged 

5-20° 

0.9 

Quasi-Newton 

yes 

5°, 10° 

0.5, 0.7 

Quasi-Newton 

yes 

O 

O 

1 

0 

0.9 

Quasi-Newton 

yes 

5° 

0.9 

Quasi-Newton 
and Least Squares 

no 

20° 

0.9 

Quasi-Newton 

yes 

20° 

0.9 

Quasi-Newton 

yes 



Table 1* (Continued). 


Wing Description 

Modeling Details 

Angle of 
Attack 

Mach 

Number 

Method of Solution 

Solution 

Converged 

63.5*^ flat delta 
wing-body com- 
bination, A = 2.0 

Without near wake, and 
with carry-over lifting 
system and the upper 
body surface separated 
by a small distance 
(0.01); no discontinu- 
ity exists between body 
and wing upper surfaces 

20‘ 

0.9 

Quasi-Newton 

yes 

78.7V63.5® flat 
double delta 
wing; A = 1.83 

Without near wake, and 
with single vortex 
system along inboard 
leading edge only 

0 

O 

CM 

0.9 

Quasi-Newton 

yes 

78.7V63.5*^ flat 
double delta 
wing; A = 1.83 

Without near wake, and 
with single vortex 
system all along lead- 
ing edge 

20“ 

0.9 

Quas i-Newton 

yes 

78.7763.5*^ flat 
double delta 
wing; A = 1.83 

Without near wake, and 
with two separate vor- 
tex systems along in- 
board and outboard lead- 
ing edges 

20“ 

0.9 

Quasi-Newton 

no 

78. 7763.5" flat 
double delta 
wing; A = 1.83 

Without near wake, and 
with two separate vor- 
tex systems along in- 
board and outboard lead- 
ing edges 

5“-20“ 

0.9 

Least Squares 

yes 



Table *1 . 


Wing Description 

78.7V63.5" flat 
double delta 
wing body combi- 
nation, A = 1.83 

78.7V63.5^ flat 
double delta 
wing-body combi- 
nation A = 1.83 


60® flat delta 
(DM-1 glider 
with thicknes s 
neglected) with 
flat leading 
edge extension, 
A = 1.8 

60® flat delta 
(DM-1 glider 
with thickness 
neglected) with 
flat leading 
edge extension, 
A = 1.8 


Modeling Details 


Without near wake or 
c ar r y-o ver lift ing 
system, and with single 
vortex all along lead- 
ing edge 


Without near wake, or 
c ar r y-ove r lift ing 
system, and with two 
separate-vortex systems 
along inboard and out- 
board leading edges 


With near wake 


With near wake 


o 


e of 


Mach 

Number 


So lut ion 

Method of Solution Converged? 


20® 0.9 Least Squares no 


20® 0.9 Least Squares no 




Table 2. 


Comparison of aerodynamic characteristics of EM-1 glider at a = 
15* and M = 0. 



c 

c 

c 

Method 

L 

D 

m 

FVS method 

0.68 

0.18 

0.010 

Data (ref. 14) 

0.55 

0.12 

-0.038 
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a, degree 


Figure 2 


C 

m 


Aerodynamic characteristics of A = 2.0 flat delta wing-body 
bination at M = 0.5. 


com- 




a, degree 


Figure 3 


C 

m 


Aerodynamic characteristics of A = 2.0 flat delta wing-body com- 
bination at M = 0.7. 







DATA (Ref. 13) 

WING-BODY COMBINATION WITHOUT CARRY OVER SYSTEM SEPARATED 
FROM UPPER BODY SURFACE BY A SMALL DISTANCE 

WING-BODY COMBINATION WITHOUT CARRY OVER SYSTEM 
WING ONLY 



UPPER SURFACE 





Q DATA (Ref. 13) 

WING-BODY COMBINATION WITH CARRY OVER SYSTEM COINCIDING 

WITH UPPER BODY SURFACE 



Figure 6. Effect of different types of modeling on spanwise pressure dis- 

distributions for A = 2.0 flat delta wing at ^ = 0.85, a = 20°, 
and M = 0.9. ® c^. » ’ 



DATA (Ref. 13) 

WING-BODY COMBINATION WITH CARRY OVER SYSTEM 
SEPARATED FROM UPPER BODY SURFACE BY A SMALL DISTANCE 

WING-BODY COMBINATION WITHOUT CARRY OVER SYSTEM 
WING ONLY 





LOWER SURFACE 


L 

Figure 7 


Spanwise pressure distributions for A = 2.0 flat delta wing-body 
combination at a - 5* and M = 0.7. 



LOWER SURFACE 


Figure 8 


Spanwise pressure distributions for A = 2.0 flat delta wing-body 
combination at a = 10® and M = 0.7. 
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Figure 17. Spanwise pressure distributions for A = 1. 83 flat double delta 
wing at $— = 0.53, a = 10°. and M = 0.9. 
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Figure 18. Spanwise pressure distributions for A = 1.83 flat double delta 
wing at ^ * 0.71, a = 10*. and M = 0.9. 
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Figure 19. Spanwise pressure distributions for A =* 1.83 flat double delta 
wing at —— - 0.89, a = 10°, and M = 0.9. 
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FVS METHOD (Ref. 1) 
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O DATA (Ref. 13) 

FVS METHOD (Ref. 1) 
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— x/c^ = 0.8'9' 

0 DATA (Ref. 13) 

FVS METHOD (Ref. 1) 
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Figure 23. Spanwise pressure distributions for A = 1.83 flat double delta 

wing at ^ = 0.89, a = 15*, and M = 0.9. 
r 
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VORTEX ON FLEE LEADING EDGE 



Figure 2$». DM-1 glider, with vertical fin removed and FLEE attached, modeled by 
FVS method. 





